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Representing an ensemble of random lasers with an ensemble of random matrices, we compute 
average number of lasing modes and its fluctuations. The regimes of weak and strong coupling of 
the passive resonator to environment are considered. In the latter case, contrary to an earlier claim 
in the literature, we do not find a power-law dependence of the average mode number on the pump 
strength. For the relative fluctuations, however, a power law can be established. It is shown that, 
due to the mode competition, the distribution of the number of excited modes over an ensemble of 
lasers is not binomial. 
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I. INTRODUCTION 

Random lasers with coherent feedback (see Ref. [1] 
for a review) are systems based on active disordered 
materials where self-sustained radiation modes can be 
formed. Another possibility involves sufficiently open 
wave-chaotic resonators filled with active media. In the 
absence of pumping, both realizations are characterized 
by short-lived strongly interacting passive modes. In or- 
der to treat such systems, the standard laser theory [2, 3], 
which assumes an almost closed resonator, needs to be 
modified. The first step in this direction was made in 
Ref. [4], where the Langevin formalism was adapted to 
the present situation. The authors described lasing-mode 
oscillations by non-Hcrmitian matrices and derived an ex- 
pression for the linewidth. In Ref. [5] a connection was 
made between the Langevin and the master equations. 
An improved treatment of non-linearities in the multi- 
mode laser theory was proposed in Ref. [6]. 

An interesting problem (also from experimentalist's 
standpoint) is the effect of mode competition on the num- 
ber of lasing modes. Equations yielding this number were 
derived in Refs. [7, 8] for a weakly open resonator. The 
authors of Ref. [8] discovered that the average number 
of modes varies (within certain limits) as a power of the 
pump strength. The theory is based on the fact that the 
passive-mode widths for a weakly open resonator follow 
a x 2 distribution. The results were tested numerically 
by generating ensembles of possible widths according to 
this distribution. In Ref. [9] the mode-number equation 
of Ref. [8] was rederived for a resonator with overlap- 
ping passive modes. In particular, Ref. [9] exploits a 
possibility to reproduce statistical properties of modes 
in open chaotic resonators with the specially chosen en- 
sembles of non-Hermitian random matrices. The average 
mode number was computed from sets of passive-mode 
widths obtained directly from these ensembles. It was 
claimed [9] that this number has a power-law dependence 
on pumping in a resonator strongly coupled to the envi- 
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ronment. The theory of Ref. [9] relies on the eigenvector 
and eigenvalue statistics of non-Hermitian random matri- 
ces, which was extensively studied in the literature [10- 
13]. However, general analytical expressions for the rel- 
evant correlations of the left and right eigenvectors are 
still unknown. The correlations were carefully studied 
numerically in Ref. [9]. 

In the present paper we employ the approach of Ref. [9] 
to study the average number of excited modes and, for 
the first time, its fluctuations. In the weak-coupling case, 
our results for the average agree with the predictions of 
Ref. [8]. For strong coupling, contrary to the results of 
Ref. [9] , we do not find a power law for the average, but do 
find it for the relative fluctuations. Relating the average 
number of lasing modes to its fluctuations, it is possible 
to extract some information about the distribution of this 
number over an ensemble of lasers. We argue that a 
binomial distribution gets distorted by the passive-mode 
overlap and the active-mode competition. 

II. THEORETICAL BACKGROUND 

We begin by recalling the derivation of Eq. (32) [8, 9], 
which yields the number of lasing modes for a given pump 
strength. This equation is analyzed numerically in the 
subsequent sections. 

A. Langevin equations 

The system under consideration comprises an open res- 
onator filled with J\f identical atoms. A random laser can 
be modeled by a resonator with an irregular shape, such 
that its eigenfunctions are chaotic. The system Hamilto- 
nian 

H sys = Hf + H a + Hf. a (1) 

is a sum of the radiation-field Hamiltonian Hf, atomic 
Hamiltonian if a , and their interaction -fff_ a . 

It is convenient to represent the system as field and 
atoms in an ideal (isolated) resonator interacting with 
environment (heat and pump reservoirs, or baths). The 
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reservoir degrees of freedom are then eliminated. The 
reservoirs acting on the field and on the atoms are as- 
sumed to be independent of each other [14]. Accordingly, 
the field Hamiltonian can be written in the form 

H f = n^2u x a\ax + AH U (2) 

A 

where a\ are the annihilation operators for the modes 
of the closed resonator with frequencies uj\ and AHf in- 
cludes the bath Hamiltonian and the resonator-bath in- 
teraction. In the case of an empty resonator with an 
opening, the role of the reservoir is played by an ex- 
ternal field having a continuous spectrum. This model, 
adopted also for the present paper, was carefully studied 
in Refs. [15, 16], where, in particular, the ways to split 
the field into the internal (resonator) and external (bath) 
parts were discussed. 

The atoms will be approximated by their two active 
levels separated by energy hv. Given the Fermionic op- 
erators ci,2 for the levels, one can define the pseudospin- 
1/2 operators a = s x — is y — c\c 2 and s = s z — 
\ (c\c2 — c\cij , i.e., Si generate a spin su(2) algebra. 

The spin operators a p and s p (p = 1, . . . , A/") of differ- 
ent atoms commute. The atomic Hamiltonian becomes 

AT 

H & = hv^s p + AH & , (3) 
P =i 

where the reservoir operators are contained in Ai7 a . 
The interaction Hamiltonian 

Jff.a = ifrj^ (gxpA^p - h - c -) ( 4 ) 

Ap 

is written in the rotating-wave and dipolc approxi- 
mations. The former neglects the terms proportional 
to a\tj p and a\<7^ p . Since v ~ uo\, such antiresonant prod- 
ucts would oscillate with a double optical frequency in the 
interaction picture. The latter approximation can be ap- 
plied since the optical wavelength is much larger than the 
atom size. Then the coupling constant (in the Gaussian 
units) 



9\ P = v^j — d 2 i • (p* x (r p ) (5) 

is expressed in terms of the dipole moment d 2 i for the 
1 — ► 2 transition between the atomic levels, as well as 
the normalized vector- valued eigenfunction cf) x (r p ) of the 
mode A at the atom position v p [17]. In chaotic res- 
onators the values of an eigenfunction at any two po- 
sitions arc uncorrelated (apart from normalization and 
boundary effects) if they are more than a wavelength 
apart [18]. Hence, the couplings g\ p can be treated as 
independent Gaussian random variables. 

In the Heiscnberg picture, an equation of motion for 
an operator O is O = j- [if, O] . For the laser operators 



oa, cr p , and s p , the equations of motion can be cast in 
the form of Langevin equations: 

a x = —icu\ci\ - ^7aa'»A' + y^flAp'<V + F\, (6) 

A' p' 

& p = - {if + 7_l) (j p + 2 ^ 9x' P ax'S p + F p CT , (7) 

A' 

s p = 7|| (S - s p ) - (ffA'p 4' °p + h - c -) + F p- ( 8 ) 

A' 

The reservoirs enter the equations via the damping (7AA' , 
7_L, 7|| ) and pumping (S) parameters and the operators of 
stochastic forces (F\, F£, Fp). The latter have zero reser- 
voir average and are <5-correlated in time. This property 
is a consequence of the Markov approximation, which re- 
quires the reservoir relaxation time to be much smaller 
than all the other time scales. Equations (6)-(8) are ap- 
propriate for chaotic resonators. They differ from the 
standard equations of the laser theory [3] in to aspects: 
the non-diagonality of 7a a' and the randomness of c/Ap- 
Equation (6) was derived in Refs. [5, 15]. The damping 
matrix 7aa' is Hcrmitian. It is strongly non-diagonal if 
the resonator modes are overlapping. The off-diagonal 
elements point to the interaction between the respective 
modes via a coupling to the continuum. Equations (7) 
and (8) follow, e.g., from the Langevin theory for three- 
level atoms [2] if the total population of the two active 
levels is kept constant. 7^ and 7|| are the polarization 
and inversion decay constants, respectively. 

In the following, we work in the classical approxima- 
tion, whereby the noise forces are neglected and the op- 
erators are treated as c-numbers (the prior notation will 
be retained). This approximation fails near the lasing 
threshold, where the average intensity is smaller than its 
quantum fluctuations. The classical version of Eqs. (6)- 
(8) becomes 

j t \a) = -iO|a)+^|<v)<V, ( 9 ) 
p' 

& p = -(w + 7_l) a p + 2 (g p \a)s p , (10) 

s p = 7|| (S - s p) - (H.9p) ct p + c.c.) . (11) 

Here, for the compactness of notation, we introduced the 
classical vectors |a) = (. . . , a\, a\ + i, . . .) T and \g p ) = 
{■■■, 9x P , 9x+i, P , ■ ■ -) T and the matrix 

(^)aa' = w a<5aa' - il\\>- (12) 

Since 12 is not Hcrmitian, it has different left and right 
eigenbases (Lk \ and \Rk), respectively. Thus, its spectral 
decomposition is of the form 

n = '$2\Rk)(to k -iK k )(L k \, (13) 
k 

where the eigenvectors are normalized in such a way 
that (L k \R k >) = 5 kk i and (R k \R k ) = 1, but, in general, 
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(R k \Rk>) ^ 8 k k>- f^fe - i^k (n k > 0) are the complex 
cigenfrequcncies of the passive open resonator. We ex- 
pand | a) = ^2 k ak\Rk), where the amplitudes 

otk = (L k \a), (14) 

satisfy the equations 

dfc = - (ififc + Kfe) a fc + ^(-LfelfiTp') V' ( 15 ) 

p' 

<t p = - (w + 7_L)cr p + 2^2(g p \R k ') a k 'S p , (16) 

fe' 

«p = 7|| (S - s p ) - (( R k> \g P ) a* k ,o- p + c.c.) (17) 
fe' 

following from Eqs. (9)-(ll). It is worth emphasiz- 
ing again that the open-resonator modes k are coupled 
through the interaction with the atoms, while the closed- 
resonator modes A, in addition, interact via the reservoir. 

B. Number of lasing modes 

We solve Eqs. (15)-(17) by treating the interaction with 
the atoms perturbatively. Namely, it will be assumed 
that the sustained field oscillations are proportional to 
the unperturbed eigenvectors \R k )- However, the oscilla- 
tion frequencies O/j « Ofe are to be determined selfcon- 
sistently. It was argued in Ref. [19], that a multimode 
solution is possible if 
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where Af2 is a typical spacing between the lasing-mode 
frequencies. This condition ensures that the population 
inversion s p is approximately constant in time [21]. An 
expression for s p can be derived as follows. First, one 
represents the polarization a p — ^ fc a p k as a sum of 
single-frequency components a pk oc exp (— iQ k t). Using 
Eq. (16) with s p = const, 



o-pk = 



2s p (g p \R k ) 
-i {Q k - v) + 7-L 



Oik 



(19) 



is expressed in terms of a k , which oscillates with the same 
frequency. Finally, these a pk are substituted to Eq. (17) 
yielding 



Y p = l + 



(20) 

— y2\(g P \R k )\ 2 £ k \a k \ 2 , (21) 

7J-7|| V 



(22) 



where the oscillating products a* k a k ,, k ^ k', were av- 
eraged out, in line with the constant-s p approximation. 
Taking into account Eqs. (19) and (20) and keeping only 



the terms oscillating with frequency Cl k in Eq. (15), we 
arrive at an equation for a k , 



a k ! 
B k 



(iQ k + Kfe - Bk) a k , 
2S 



(Q k - v) + 7_l 



(ifeiE 



o)(g P 



(23) 

«• (24) 



Equations (23) for N lasing modes are equivalent to a 
system of 2 N real equations [9] , 



flk = fl k — lmB k , 
ReB k = Kk, 



(25) 
(26) 



from which the frequencies Sl k and the intensities I k — 
\a k \ 2 of the lasing modes can be determined. Equa- 
tions (25) and (26) are valid only for such k that I k > 0. 

Further progress can be made if Y~ 1 in Eq. (24) is 
expanded up to the linear terms in I k . This procedure 
presumes that the laser is operating not far from the 
threshold. When the atoms are distributed uniformly 
over the resonator and their density is sufficiently large, 
the p sum becomes a volume integral. Then the summa- 
tions entering Eq. (24) are computed as follows: 



J2\9p)(9p\^^9 2 / d 3 r\d>*(r))(<j>*(r)\ = Mg 2 , (27) 
p J 

^2(Lk\g P )(g P \Rk)\(g P \Rk')\ 2 
p 

^Mg A V j d?r{L k \<i>*{T)){4>*{v)\Rk)\{cl>*{v)\Ry)\ 2 
^Ng\l + 28 kk ,). (28) 

where \c/>(r)) = (. . . , (f>\{r), <t>x+i (r), . . .) T , 



2-kv 



121 1, 



(29) 



and V is the resonator volume. Above we assumed that 
all the modes A are polarized along d2i and u\ « v. 
The approximate equality in Eq. (28) results from treat- 
ing the wavefunctions (L k \(j)*(r)) and ((p*(r)\R k ) in a 
wave-chaotic resonator as Gaussian random variables, 
restricted only by normalization. Performing random- 
matrix simulations, this property was shown to hold in 
the relevant range of the mode widths [9]. Using Eqs. (27) 
and (28), we arrive at 



2SA/y / n k 

B k = £fc 1 + 1- 



7± 
1 



4ff 2 



7±7|| k , 
The linear gain Go£ k , where 



7-L 

I k 'C k '{l + 25 kkl ) 



+ Q(I 2 ). 



7-L 



(30) 



(31) 



is obtained from Re by setting = 0. 

To determine the number of lasing modes, we approxi- 
mate fifc with Ofc, substitute Bk (30) in Eq. (26), divide it 
by Ck, sum over k, and find J2k' Ik'£-k>- Then J2k' h'C-k 1 
can be used in Eq. (26) to express Ik , which is required to 
be positive for all lasing modes. This condition yields [8] 



(32) 




where the modes are ordered in such a way that K\/ C\ 



< 



< 



The largest N satisfying this inequality is 



the number of lasing modes, 
largest iVo such that 



On the other hand, the 



£n 



<G 



(33) 



holds, is the number of modes that would lase in the 
absence of mode competition. These are the modes 
for which the linear gain exceeds the losses. Clearly, 
N > N. 




FIG. 1: (Color online) Distribution p(fi) of the real parts of 
the eigenvalues of random matrices. Numerical simulations 
were performed for an ensemble of 500 matrices of size L = 
200 coupled to M = 0, 1, 5, and 10 open channels with the 
coupling 7 = 1. Analytical expression (35) is plotted for 
comparison. 



III. AVERAGE NUMBER OF LASING MODES 
AND ITS FLUCTUATIONS 

A. Random- matrix model 

We investigated the mode statistics resulting from 
Eqs. (32) and (33) . Ensembles of open chaotic resonators 
were modeled by randomly generated 11011- Herrnitian ma- 
trices ft = ilt — (12). In the basis of the modes A, f2 has 
a diagonal real part. Its imaginary part is of the form [5] 

-^=--kWW\ (34) 

where (W) Xm (A = 1, L, m = 1, M < L) de- 
scribes interaction of the resonator mode A with the mth 
channel of the reservoir. Each channel includes a contin- 
uum of frequencies. However, W is frequency indepen- 
dent in the Markov approximation. Clearly, the matrix 
7 has at most M nonzero eigenvalues 7^, i = 1, . . . , M (it 
is possible to find N — M linearly independent vectors 
orthogonal to the rows of W^). We will consider the case 
of equivalent channels, when all 7; =7. 

A matrix Cl is most easily constructed in the basis 
where 7 is diagonal: we fix 7 and choose u> from a Gaus- 
sian orthogonal ensemble [12]. Without loss of general- 
ity, the diagonal and off-diagonal elements of to are taken 
from a normal distribution with zero mean and the vari- 
ance of 2/L and 1/L, respectively. In the limit L ^> 1, 
£lk (13) (the real parts of the eigenvalues of Ci) are dis- 
tributed according to the Wigner semicircle law 



p{Q) = -yjl-—, -2<ft<2, (35) 

where p(fl) is normalized to unity. Numerical simulations 
(Fig. 1) show a reasonable agreement with this equation. 



The strength of coupling to continuum is characterized by 
a parameter 2irp(Q)/ (7 + 7 _1 ) [12]. Thus, it is sufficient 
to consider 7 £ [0, 1], whereby 7 = (7 = 1) corresponds 
to the vanishing (strongest) coupling. Importantly, even 
within one matrix Q, the effective coupling depends on 
the spectral region according to p(f2). 

B. Results and discussion 

In the following figures, we present results of numerical 
simulations for the average number of lasing modes (N) 
and its standard deviation a^. The respective quanti- 
ties in the absence of the mode competition, (No) and 
(TAr , were calculated as well. The averages were per- 
formed over ensembles of random matrices. As was ex- 
plained earlier, the effective coupling to continuum de- 
pends on p(Q). Therefore, for each matrix, of all eigen- 
values, only Lo ~ 0.36L eigenvalues closest to the top of 
the Wigner semicircle were used in Eqs. (32) and (33). 
Within this spectral region, p(O) varies by about 4%. In 
order to reduce the number of parameters, we assumed 
that 7^ is sufficiently large, and set Ck = 1. The pump- 
ing was measured in units of its threshold value 5thD 
which was determined numerically from the threshold 
condition (N) = 1. An estimate yields S/Sthr ~ Go/no, 
where ko ~ -fM /L is a typical loss. 

First, we discuss a weakly open resonator, kq -C A57 w 
ir/L, where Aft is the mean nearest-neighbor spacing be- 
tween fife. A theory for the averages (N) and (iVo) in this 
regime was proposed in Ref . [8] . Central to the argument 
is the analytical expression for the distribution of kj, (13), 

, , (M/2) M / 2 m, / M \ 
P(V)= r(M/2) W - exp(- T yj, y = k/k, 

(36) 
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FIG. 2: (Color online) Average number of lasing modes (N), 
its standard deviation <jjv, and relative standard devia- 
tion <jn/{N), as a function of pumping. The same in the 
absence of mode competition ({No), <jn , ctn /(No)). The 
thin solid lines display analytical results of Eqs. (37) and (38) 
and power-law asymptotics. Numerical simulations were 
performed for an ensemble of 75 random matrices of size 
L — 1000 coupled to M = 5 open channels with the cou- 
pling 7 = 0.1. Lo = 360 eigenvalues per matrix were taken 
into account. The analytical functions were computed with 
the parameters k = 5.0 x 10 -4 and Go |s=s thI = 3.3 x 10~ 5 , 
which were extracted from the ensemble. 



where k is the average of k and T(z) is the gamma func- 
tion [20]. P(y) is a \ 2 distribution with M degrees of 
freedom and the average y = 1. The average number of 
modes in the absence of competition can be estimated 
from 



(N ) = L Q 



Gq/k 



dyP(y)=L Q 



7(f. 



M Go ' 
2 re , 



r(M/2) 



(37) 



where j(z,x) = J Q X dtt z ~ 1 e~ t is the incomplete gamma 
function [20]. The numerical results in Fig. 2 show a 
good agreement with this prediction. In the weak-pump 
regime Go/k <C 1, there is a power law (No) oc S M / 2 [8]. 
With the increased pumping, the saturation (No) — > Lo 
sets in. Clearly, this form of saturation is an artifact 
of our model. Normally, the number of potential lasing 
modes would be limited by the Lorentzians Ck- Never- 
theless, for large Lo, the present model is appropriate 
below the saturation. 

In order to determine <t/v , it is the simplest to as- 
sume that No — 1 , ■ • • , Lo is distributed according to a 
binomial distribution P p (No\Lq), where p = (Nq)/Lq is 
the probability for a given mode to lase. The standard 
deviation for this distribution is known to be 



(No) 1- 



(N ) 
Lo 



(38) 



where (No) can be substituted from Eq. (37). A com- 
parison with the numerical curve in Fig. 2 supports our 



FIG. 3: (Color online) Same as in Fig. 2, but with the coupling 
7 = 1.0. The thin solid line displays a power law with the 
exponent —0.32. The ensemble is characterized by 7t = 7.4 x 
10" 3 and Go !s=s thr = 1.1 x 10" 4 . 



assumption. When (Nq) <C L , we find a power law 

In the presence of mode competition, it is interest- 
ing to look at the case when (N) is far below the sat- 
uration, but, still, sufficiently large. The former condi- 
tion yields (N) oc k m / 2 , while the latter ensures that 
the terms of order Nkn ps (N)k dominate the left-hand 
side of Eq. (32). Combination of the two estimates gives 

M 

(N) oc [8]. We checked numerically (via ctat) that 

the N distribution is not binomial. Nevertheless, the 

I M 

power law oc J (TV) oc S 2 <- M +v remains valid (Fig. 2). 

Next, we consider the problem of a strong coupling to 
the bath (kq ~ Af2) modeled here by a random-matrix 
ensemble with 7=1. The distribution P(y) in this case 
is no longer given by Eq. (36), but rather has a power- 
law tail P(y) oc y~ 2 , y ^> 1 [12]; its full analytical ex- 
pression is very complicated. Numerical results for the 
average number of lasing modes and its fluctuations were 
obtained for M = 1,3, 5, 7, and 10. The data for M = 5 
are presented in Fig. 3. The S dependencies of (N) and 
(No) here are similar to those found in Ref. [9]. How- 
ever, except for M = 1, we can confirm for these quan- 
tities neither a power-law behavior, in general, nor the 
powers 1/3, respectively, 1/2, in particular. 

An analysis of fluctuations shows that the distribu- 
tions of N and A^o are non-binomial. While the standard 
deviation er/v saturates at large S, the relative deviation 
&n I (N) exhibits a power law with the exponent 0.2-^0.4, 
depending on M. Examining the numerical data, we dis- 
covered an interesting relation between the relative fluc- 
tuations of No and the average of N: 



a N 

(No) 



1 

(AO 



1 

T : 



s 

Sthr 



> 1. 



(39) 



As demonstrated in Fig. 4, this property is well satisfied 
for M^l. Unfortunately, an explanation of this result 
is still lacking. 
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FIG. 4: (Color online) Comparison of (TV) -1 - L^ 1 (thick 
lines) and ctn /{No) (thin lines) plotted as functions of pump- 
ing. The ensembles of random matrices as in Fig. 2, but with 
the coupling 7 = 1.0 and with the number of open channels M 
labeled in the figure. 



IV. CONCLUSIONS 

By modeling ensembles of open chaotic resonators with 
ensembles of random matrices, we studied average num- 
ber of lasing modes and its fluctuations. To highlight the 
effect of mode competition and allow for a better com- 
parison with earlier work, the linear approximation to 
the lasing equations was considered as well. 

In the case of a weakly open resonator, the average 
number of modes is proportional to a power of pump 
strength (within a certain pumping range). This result 



agrees with the analytical prediction in Ref. [8]. The 
standard deviation changes as a square root of the aver- 
age. In the absence of mode competition, the number of 
modes follows a binomial distribution, if the total num- 
ber of eigenstates available for lasing is finite. The dis- 
tribution becomes non-binomial in the presence of mode 
competition. 

For a resonator strongly coupled to the environment, 
we could not establish any power-law dependence of the 
average mode number on pumping. (The case of one 
open channel makes an exception.) This evidence stays 
in contradiction to the conclusions of Ref. [9]. On the 
other hand, we find a power-law behavior for the relative 
fluctuations. A curious relation (39) between the aver- 
age with and the relative fluctuation without the mode 
competition requires further investigation. 

As a possible extension of this work, it would be in- 
teresting to relax some of the assumptions made. For 
example, one can consider an effect of the Lorentzian 
line shape (here approximated as rectangular). A more 
challenging task is to avoid the near-threshold expansion 
of the intensity-dependent denominators in Eq. (24). 

Acknowledgments 

The author is grateful to Fritz Haake for introducing 
him to the research area of random lasers and continuing 
support of the project. Hui Cao, Gregor Hackenbroich, 
Dmitry Savin, and Hans-Jiirgen Sommers are acknowl- 
edged for helpful discussions. This work was financially 
supported by the Deutsche Forschungsgcmcinschaft via 
the SFB "Transregio 12." 



[1] H. Cao, Waves Random Media 13, Rl (2003). 

[2] M. Sargent III, M. O. Scully, and W. E. Lamb, Jr., Laser 
Physics (Addison- Wesley Publ. Co., Reading, 1974). 

[3] H. Haken, Light, vol. 2 (North-Holland Publ. Co., Ams- 
terdam, 1985). 

[4] G. Hackenbroich, C. Viviescas, and F. Haake, Phys. Rev. 

Lett. 89, 083902 (2002). 
[5] G. Hackenbroich, C. Viviescas, and F. Haake, Phys. Rev. 

A 68, 063805 (2003). 
[6] H. Tiireci, A. D. Stone, and B. Collier, e-print 

arXiv:cond-mat/0605673, 2006. 
[7] H. Haken and H. Sauermann, Z. Physik 173, 261 (1963). 
[8] T. S. Misirpashaev and C. W. J. Beenakker, Phys. Rev. 

A 57, 2041 (1998). 
[9] G. Hackenbroich, J. Phys. A: Math. Gen. 38, 10537 

(2005). 

[10] Y. V. Fyodorov and H.-J. Sommers, J. Math. Phys. 38, 
1918 (1997). 

[11] Y. V. Fyodorov and B. A. Khoruzhenko, Phys. Rev. Lett. 



83, 65 (1999). 

[12] H.-J. Sommers, Y. V. Fyodorov, and M. Titov, J. Phys. 

A: Math. Gen. 32, L77 (1999). 
[13] Y. V. Fyodorov and H.-J. Sommers, J. Phys. A: Math. 

Gen. 36, 3303 (2003). 
[14] H. Haken, Laser Theory (Springer- Verlag, Berlin, 1984). 
[15] C. Viviescas and G. Hackenbroich, Phys. Rev. A 67, 

013805 (2003). 

[16] C. Viviescas and G. Hackenbroich, J. Opt. B: Quantum 

Semiclass. Opt. 6, 211 (2004). 
[17] D. F. Walls and G. J. Milburn, Quantum Optics 

(Springer- Verlag, Berlin, 1995). 
[18] M. V. Berry, J. Phys. A: Math. Gen. 10, 2083 (1977). 
[19] H. Fu and H. Haken, Phys. Rev. A 43, 2446 (1991). 
[20] M. Abramowitz and I. A. Stegun, eds., Handbook of 

Mathematical Functions (Dover, New York, 1972). 
[21] The same is tacitly assumed in Ref. [9] 



